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Abstract: The black-hole black-string system is known to exhibit critical dimensions and 
therefore it is interesting to vary the spacetime dimension D, treating it as a parameter of 
the system. We derive the large D asymptotics of the critical, i.e. marginally stable, string 
following an earlier numerical analysis. For a background with an arbitrary compactifica- 
tion manifold we give an expression for the critical mass of a corresponding black brane. 
This expression is completely explicit for T n , the n dimensional torus of an arbitrary shape. 
An indication is given that by employing a higher dimensional torus, rather than a single 
compact dimension, the total critical dimension above which the nature of the black-brane 
black-hole phase transition changes from sudden to smooth could be as low as D < 11. 
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1. Introduction 

In the presence of extra compact dimensions, there exist several phases of black objects 
depending on the relative size of the object and the relevant length scales in the compactifi- 
cation manifold. More precisely, by "black objects" we mean massive, non-rotating (static) 
solutions of General Relativity with no extra matter 1 . 

So far most research on this problem concentrated on the background IR" 0-1 x S 1 , 
where S 1 is a single compact dimension and often D has taken the values 5, 6. However, we 
know of two instances of a critical dimension in this system, where the qualitative behavior 
changes. Therefore, we consider spacetimes of total dimension D out of which d < D 
dimensions are extended, and we vary both D and d as parameters of the system. 

The evidence is as follows. In the first case ^ evidence was given that the behavior of 
the system near the conjectured point of merger between the black hole and black string 
branches depends on a critical dimension -^merger = 10, such that for D < ^merger there 
are local tachyonic modes around the tip of the cone (conjectured to be the local geometry 
close to the thin "waist" of the string) which are absent for D > L>merger- The second 
case relates to a different point in the phase diagram - the Gregory-Laflamme point (GL) 
where the string is marginally tachyonic [|[. It was recently found || that the order of 
transition has a critical dimension 2 Dq L ="13.5", such that for D < Dq L the transition is 
first order, and otherwise it is second order. Moreover, there is another example, closely 

1 We assume these are non-singular outside their event horizon. 

2 Of course dimensions are integral, and the notation means only that the change in the order happens 
between D — 13 and D = 14. 
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related to the first, of a critical dimension in a different system - the Belinskii-Khalatnikov- 
Lifshitz (BKL) analysis of the approach to a space-like singularity, where there is a critical 
dimension D* BKL = 10, such that for D > D* BKL the system becomes non-chaotic (see the 
review Q| and references therein). This last example is related to our system not only in 
that it displays a critical dimension, but also in that the equations are similar to those of 
the merger transition - the BKL equations for the metric components as time approaches 
the singularity are similar to those for the cone perturbations near its tip. 

A description of relevant works follows. Gregory and Laflamme |j| discovered the 
critical mass for a uniform string, [iql, and determined it (numerically) for 5 < D < 10. 
Actually, as was realized by Reall [||, in order to determine the critical mass it is enough 
to know the eigenvalue for the negative mode of the Euclidean Schwarzschild geometry 
discovered by Gross Perry and Yaffe (GPY) ^ and evaluated there in Ad. The analysis of 
@ was generalized to d dimensions by Prestidge [|j (see also Q for a related discussion of 



black holes in a box). Gubser [10] set a method of perturbation analysis around the 5D 
GL critical string whose first step is to determine fiGL by a somewhat different way than 
using 0. That method was perfected by Wiseman |0, |12| who also applied it in 6D. Its 



generalization to an arbitrary dimension was made in [|| , (see appendix |A| for a summary) . 
In fact, by now there is a considerably growing literature on this black-hole black-string 
system. See [12-17] which study non-uniform strings and [18-26] which concentrate on 
black holes in the system. 

In this paper we obtain some results on the dimension dependence of the GL effect 
and in particular of the critical mass. Recall that the dimensionless mass of a black 
string is defined by n := GdM/ L D ~ 3 where Gd is the D dimensional Newton constant 
and L is the asymptotic length of the compact circle. In [|| the critical mass of a black 
string on R- 0-1 x S 1 was obtained numerically for a range of -D's, up to D = 50, and 
a "phenomenological" exponential behavior ugl oc j d was observed where 7 « .68. In 
section |2| we derive the large D asymptotics of the critical mass and determine that 7 = 
y' 'e/2ir in agreement with that "phenomenology" . It is interesting to compare this behavior 
with that of the "equal areas for equal mass" estimator for the critical string 3 , ^5: both 
expressions are exponential in D but with a different base since for D > 1, ^5 a ^y s D 
where 75 = e -1 / 2 ~ 0.606. 

Moreover, the behavior of the critical dimensionless mass turns out to be equivalent 
to a simple behavior of the critical Schwarzschild radius: ro ~ (L /2-n)^/d, showing that 
for large d the critical string is actually "fat" and as we show in section ||| most probably 
it cannot decay into a black hole. 

In order to derive the asymptotics we use the equation for the negative mode of the 
d dimensional Euclidean Schwarzschild geometry which is a single, second order, ordinary 
differential equation (ODE). In the large d limit the equation becomes the radial part of 
the flat space Laplacian with all the non-trivial information being carried by the boundary 
condition close to the horizon. 



A This is defined to be the mass for which the areas of the uniform black string and a localized, naively 
undeformed black hole are equal. 
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In section ||] we proceed to vary the dimension of the compactification manifold and 
we express the critical mass of a uniform brane in this general background in terms of 
the lowest eigenvalue of the Laplacian operator on that manifold. In particular, when 
the compactification manifold is the flat n-dimensional torus of an arbitrary shape the 
expression can be made explicit in terms of the shortest vectors in the reciprocal lattice. 
We note that there are interesting cases when several modes turn tachyonic simultaneously. 

Finally in section || we raise the question whether it is possible to see the second order 
behavior for dimensions smaller than 14 using T ra compactifications. It would be especially 
interesting if it happened for dimensions D < 11 where M-theory/string theory is believed 
to be a consistent theory of quantum gravity. As in || we use as an estimator for the 
critical dimension the intersection between hgl(D) and ns{D). We find that the critical 
dimension can indeed be reduced and for 3 < n < 6 it is estimated to be around D = 10 
(see figure . 



2. Large dimension asymptotics for ugl 

Consider a D dimensional background with d < D extended dimensions, and the other 
D — d being compact. In this section 4 d = D — 1 and the compactification manifold is S 1 
of period L. 

A perturbation theory around the static uniform black string in D dimensions was 
recently employed in |3| to construct the non-uniform string branch emerging from the 
critical GL point. Among other results it was observed that the critical mass follows 
essentially an exponential law function of D 

VOL oc J D , (2.1) 

where numeric analysis yields 7 ~ 0.686, the prefactor is approximately 0.47, and our 
definition of the dimensionless mass is as usual 

H := G D M/L D ' 3 , (2.2) 

where Gr> is the D dimensional Newton constant 5 . 

This critical mass was compared with us, the commonly used "equal areas for equal 
mass" estimator. More precisely it is an estimator for the point of the first order transition, 
and it is only an estimate since in the absence of exact solutions for caged black holes 6 one 
approximates the mass-area relation with that of a small spherical black hole. One finds 

1 t] D °3 3 (j-3)(^)(M 
MS ~ 16tt Q^-i (D - 2 )P-2)(fl-4) ' 

where D _i := D 7r D / 2 / (D/2)\ is the area of the unit S D 1 sphere. 



4 This assumption will be relaxed in later sections. 
5 In this paper we work in units such that Gd = 1. 

6 By a "caged black hole" we mean a black hole in a compactified spacetime. 
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Performing a series expansion for D — > oo we obtain 

"log(D) 



log(/i S ) = ~D ( 1 



(2.4) 



Hence, for large D, fis also exhibits an exponential scaling, though with a different ex- 
ponent, 7s = l/\fe 0.606. Since 7s < 7, fis becomes much smaller than fiQL as D 
increases. Moreover, there is a dimension, D ~ 12.5 for which us intersects This is 

an indication that for higher dimensions the black string becomes unstable before the black 
hole phase can have superior entropy. This, in turn, can be regarded as a hint that above 
that dimension the unstable GL-string decays to a state different from the black hole, and 
it is plausible to expect it to decay into a slightly non-uniform string. And indeed, exact 
calculations in higher orders of perturbation around the GL-point reveal an existence of a 
critical dimension (D* = 13.5) where the order of the phase transition between the uniform 
and the non- uniform black strings solutions changes from the first to second order Q. 

It is desirable to have an analytic formula for the numerical constant 7 in ( |2.1| ). Un- 
fortunately, the equations that were solved in Q in order to obtain the scaling ( |2.1D are 
two coupled second order, rather complicated, ODEs that required a numerical treatment. 
In the large D limit these equations simplify, but nevertheless we did not find an analytic 
solution. Hence, motivated by Reall's realization ||, we turned to look at the negative 
mode of the Euclidian Schwarzschild geometry, which yields a single ODE which we were 
able to solve analytically in the large D limit. 

2.1 The negative mode of the ti-dimensional Schwarzschild 

We aim to identify the negative modes of perturbed Euclidean Schwarzschild- Tangherlini 



metric [27] in d (extended) dimensions. The line element for that background is 



ds 2 = +f(r)dt 2 + /(r)- 1 ^ 2 + r 2 dn d _ 2 , 

ttr) = l-(^Y~ 3 , (2.5) 

t is periodic with period 4nro/(d — 3), and in this section we work in units where 7*0 = 1. 

The (gauge invariant) spectrum of the perturbations (f) a b in the transverse-traceless 
gauge is obtained by solving the eigenvalue equations 

AzAfe = A0 afe , (2.6) 
where in full index notation the Euclidean Lichnerowicz operator is A^ ab crf = — 8^8f □ — 

For this spacetime the negative eigenvalue perturbation is static and spherically sym- 
metric 0, § . Let us write the ansatz for such a perturbation as 

4> a = diag-j>(r), x(r), n(r), K(r)} . (2.7) 

* v ' 

d—2 terms 

The traceless condition sets 

< r ) = -J^{x(r)+Hr)} (2-8) 
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while transversality, 



ip(r) 



0, implies 
2rf 



r f + 2(d-l)f 
X (r) + t, 7T* X\r) 



rf'-ir y ' rf> - 2f 

Hence, the equations ( [2.6| ) reduce to a single linear second order ODE for x( r ) 
"2r 2 (ff»-f>2)-r(d-2)ff' + 2df- 



(2- 



/X"(r) + 



+ 



r(r/'-2/) 
r 2 / / / // + r [2( d _i) // //_( d + 2 )/«] 



4//' 



r(rf'-2f) 



X (r) = -k 2 X (r) (2.10) 



where for convenience we define A = —k 2 , such that for the negative eigenvalues k is 
real. The 4d version of this equation was analyzed in j?]], while in || (see also @) the 
generalization to arbitrary d was considered. 

This equation has three singular points in the region outside the horizon r > tq: at 
the horizon r = tq itself, at infinity and at r = r s which is the solution to rf — 2f = 0. r s 
is given by 

r s d ~" = ^ (2.11) 



and we note that it is exactly the critical radius for trapping light, namely, a light ray 
originating from infinity will fall into the black hole if (and only if) it crosses r s . Regularity 
of x{ r ) a t tli e horizon (a regular singular point) and at infinity (an irregular singular point) 
determines the solution. At the horizon, where f(ro) = the characteristic exponents are 
— 1, namely x ~ ( r — r o) a with a = —1, 0, and so regularity of (2.10) imposes 

k 2 



d 



2(d-3) 



(2.12) 



At r = r s , on the other hand, the characteristic exponents are a = 0, 3 and both solutions 
are finite (though with possible log's). Looking at a generic solution with x ~ (r — 



find 



X'(r s )/x(r s 



1 



._[r/' + 2(d-l)/] i 



d_ 



(2.13) 



It is important to realize that this relation is not a boundary condition, but rather a 
consequence of the equation itself. We observe that since the characteristic exponents are 
0, 3, regularity at r s can be obtained by extending the equation ( |2.10| ) to become fourth 
order (such that the exponents a = 1,2 are solutions as well), which is what happens in 
Gubser's gauge (see the Appendix). 

2.2 Large d limit 

In general, the solution to the equation ( 2.1C| ) must be obtained numerically [7-9]. However, 
let us consider this equation in the limit d — > oo. In this limit f(r) — > 1 and f'(r), f"(r) — ► 
for r > r s > 1 and hence the eigenvalues equation ( 2.10| ) simplifies to 



X"(r) + - X '(r) - k 2 X {r) = , 



(2.14) 
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which is the flat space Laplacian (in a d + 2 dimensional spacetime) . In this approximation 
r s is not singular anymore and hence ( 2.1 3| ) does not hold automatically. Instead we impose 
it as a boundary condition 

X'(r s )/ X (r s ) = -- • (2.15) 
r s 

A more rigorous treatment would derive this effective boundary condition by solving the 
original equation in the vicinity of r = 1, r s and matching it in the large d limit with the 
solutions of ( [2.14D . 

We shall now show that the solution to ( |2.14 ) subject to the boundary condition ( 2.15| ) 
and regularity at infinity can be found analytically. 

Changing variables in ( [2. 14 ) by x( r ) = u(r)r~ v with v = (d — l)/2 we obtain for u(r) 

r 2 u"(r) + ru(r) - (v 2 + k 2 r 2 )u(r) = (2.16) 



which is just a modified Bessel equation whose solution is |28|] 

u{r) = d K v {kr) + C 2 I„{kr) (2.17) 

where l v and K v are the modified Bessel functions of the first and the second kind re- 
spectively pq[ . The asymptotic boundary condition forbids a growing exponential so- 
lution. Hence 6*2 = and the selected solution is K u , whose asymptotic behavior is 
K v {z) ~ e~ z /\/2irz. The boundary condition (2.15) dictates 

xK' u (x) + {d-v)K v {x) = 

x := kr s . (2.18) 

Using the recurrence relations for K v one arrives to a simpler algebraic equation 

K v {x) -xK v ^{x) = , (2.19) 

that determines the eigenvalue k as a function of v. We note that in order to arrive at this 
equation one relies on the cancellation d—2v = 1, which is a delicate feature of this process. 
The leading behavior of k = k(d) can be gotten by making the self-consistent assumption 
k d, in which case one may use the leading behavior K v {x) ~ T(i/) (x/2)~ u /2 valid for 
x -C v. Substituting into ( 2. 1S| ) we get 

k = —Vd-3~Vd , (2.20) 

r s 

confirming the assumption k d, and at the same time showing that the critical string is 
also surprisingly "fat" vq/L = k/(2ir) S> 1 at large d. 

Eq. ( [2.19D may be solved as a series in 1/d although it is probably not justified since it 
represents only the leading d limit of the original equations where various 1/d corrections 
were already neglected. Still we record that if one uses more orders in the Taylor series of 
K v [x) the solution of ( ^l9| ) can be expanded as x = Vd - 1/Vd + 0(l/d 3 / 2 ) for d > 1. 
This formula yields less than 0.2% discrepancy with the numerical solution of ( |2,19| ) for 
d = 20 and the discrepancy is even smaller for larger d. Moreover, this estimate differs by 
less than 4% from the full solution of eq. ( |2.10 ) for d = 50 and by 2.5% for d = 100 (see 
Table |l] in appendix [A].) This is an additional confirmation of the relationship between the 
GL tachyon and the GPY negative mode. 
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2.3 The critical mass 

We can get now some insight into the mass scaling ( |2,1| ) of the marginally stable black 
strings. From the definition of \i Q2.2| ) with L = 2-K/kcL we get the critical mass for the 
GL-instability 

(d-2)n d _ 2 G D d ^fk GL \ D ~ 3 (d-m^i (r k GL \ d - 3 

(recall that Gd/G^ is the volume of the internal manifold which is L in this section). 
Expanding \og(nGL) m series for d — > oo and keeping the leading terms we get 



log (vGL(k GL )) = d log ( ) ( 1 + O 



\ log te^ + log(d) 



(2.22) 



Since according to ( |2.20| ) k G ^ — > \fd in this limit, we obtain from the above expansion 

log(/i GX ) ~ d log {J^j (2-23) 
Hence we learn that for large D = d + 1 the critical mass should scale as \i oc 7^ with 

7 = Jl. 0.658 (2.24) 

V 27T 

The numerical value 7 = 0.686 differs from this 7 by about 4%. Note, however, that the 
numerical estimate was obtained for relatively small Z)'s Q, which should account for this 
deviation. 

3. Torus compactification 

In this section we determine the critical mass for GL instability of black branes on the 
D-dimensional background M. d x T n , namely d extended dimensions and n compactified on 
a torus, and D = d + n. The torus is completely general, so its period vectors may be of 
any size and with arbitrary angles between each other. 

The relevant black branes are the uniform n-branes i3 n =SchwxT n , namely those 
which do not depend on the internal coordinates z = z l , i = 1, . . . ,n. Their metric is 



given by a sum of the d dimensional Schwarzschild metric (2.5) and the flat torus metric 
ds^n = dz % dz 1 . In addition, the torus satisfies the periodicity boundary conditions: it is 
defined by a lattice which may be given by n period vectors e% such that any function Y(z) 
must satisfy Y(z + e) = Y{z) where e belongs to the lattice: e = Y17=i m * ^ W1 ^h m ! being 
arbitrary integers. 

The critical mass is defined such that a perturbation mode h^ iu (r,z) turns marginally 
tachyonic 

A L h^{r,z) = 0, (3.1) 

where A/, is the Lichnerowicz operator in the total space B n . Separating the compact 
variables from the extended ones, and restricting to perturbations which are scalar with 
respect to T n we have 

h^(r,z) = h flu (r)Y(z) . (3.2) 
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The Lichnerowicz operator in the total space decomposes into 



A L h^(r, z) = Y(z) A L h^(r) + h^(r) (-AY(z)) , (3.3) 

where A is the Laplacian operator in the compact space T n and A^/i /iI/ (r) describes the 
Lichnerowicz operator on the d dimensional Schwarzschild geometry. Denoting the eigen- 
values by Ag chw , At™ 

A L V( r ) = A Schw V( r ) 
-AY(z) = X T nY(z) , (3.4) 



we have that the total eigenvalue (which should vanish by 3A) can be written as the sum 



= A Schw + A T' 1 • ( 3 - 5 ) 

In order to find a zero eigenvalue one of the eigenvalues should be negative. Actually 
there is only one such mode - it is the GPY mode of the Schwarzschild geometry ^ . For 
this reason it is enough to consider modes which are tensor on the Schwarzschild geometry 
and scalar on the torus. So we take 

k 2 

A Schw = —Ti ( 3 - 6 ) 



where kj 2 are the eigenvalues defined in ( 2.1C| ), computed numerically in || and whose 



large d asymptotics is Q2.2C ). Substituting back into (3.5) we find that the critical rg below 



which the black brane destabilizes is 

r 2 = ^~ (3.7) 
A min 

where A m j n is the minimal (non-zero) eigenvalue of the Laplacian on T n . Actually the 
argument is more general and applies to any stable compactifying manifold. In the case of 
T n this eigenvalue may be described explicitly 

A T - = \k\ 2 (3.8) 

where k is in the reciprocal lattice, namely k ■ e = 2n j, j £ 'K for any lattice vector e. 
Finding the minimal eigenvalue of the Laplacian reduces to finding the shortest vector in 
the reciprocal lattice of T". 



The qualitative options 

There are interesting (even if degenerate) cases where there are several vectors in the 
reciprocal lattice which lie closest to the origin, for example, a cubic lattice and a triangular 
one. Clearly that would mean that several modes would turn marginally tachyonic at the 
same brane mass 

# (marginal tachyons) = 
= # (reciprocal lattice vector closest to the origin) := : |fc| 2 = l^l™^}/^ , (3.9) 
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Figure 1: The shaded area designate the modular domain of T 2 (the sides I and I' are identified). 
The two tachyon instability develops for torii inhabiting the arc r = exp(i8) with ir/3 <0< 2ir/3. 

where in the last expression k is in the reciprocal lattice and the division by 2 accounts for 
the trivial degeneracy of k and —k. 

An example: T 2 

For the example of T 2 we can find the locus of these degenerate torii explicitly. They 
are precisely the torii such that the basis vectors e\ , e% satisfy \e\\ = | e*2 1 and the angle 
between them, 9, is in the range tt/3 < < tt/2. In terms of the modular parameter r that 
means that it lies on the boundary of the modular domain on the arc r = exp(i8) with 
vr/3 < 6 < 2tt/3 see figure ([[]). 

4. Indication for critical dimensions with T n compactification 

In the black-string case there is a critical dimension D* = 13.5 for the phase transition 
between the uniform and the non- uniform states. Namely, for a compactification on T , 
the transition is of first order below D* , while above D* it is continuous. In || it was 
shown that D* is well estimated by D* defined by the intersection [1gl{D*) = pis{D*). Let 
us generalize and estimate such a critical dimension for the T n compacification. At the end 
of this section we mention other estimates which inquire whether at given [i a black-hole 
can "fit into" the compact dimension (as stressed recently in [p9f] ). 

First, considering a single tachyon instability, when one of the torus dimensions is 
much larger than all the rest, then the situation is exactly analogous to that of the black 
string transition. Hence, in this case the critical dimension is still 13.5. 
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Figure 2: The dimension D* for which hgl{D* ,n) = ^s(D*,n) as a function of the dimension n 
of the square torus T ra . The closest integer dimension above D* estimates a critical dimension for 
a change of order in the black brane phase transition. For n = 1 this estimate is known to be short 
by about 1 from the actual D*: DJ ~ 13.5 ~ DJ + 1. For 3 < n < 6 the estimate is around 

J [71—1) [71=1) — — 

10D, making it plausible that the actual critical dimension may be as low as D* < 11, where a 
consistent theory of quantum gravity is believed to exist. 



Another important example is that of a square n-torus where, as we found in previous 
section, n tachyons appear simultaneously below the critical mass. In this case one expects 
the unstable modes to deform the horizon along all internal directions. We would like 
to know whether the transition to the non-uniform state is smooth or not. Therefore, 
generalizing equation ( |2.3j ) from the string case, we define an "equal area for equal mass" 
estimator 



1 

Ms 



16vr 



(d _ 2) (i?-3)(d-2) 1 1 " 

n£z! (d - 2p- 2 )( d - 3 ) 



(4.1) 



The estimate for a critical dimension is obtained by such a D* that ^s(D*,n) = 
Hgl{D* ,n) = /jgl(D* — n). In figure [| we plot this D* as a function of the dimension of 
the internal torus space, n. Above the curve the black hole state is not entropically favorable 
at the critical mass. This indicates probably a smooth transition to a non-uniform brane 
emerging from the GL-point. Interestingly, the value of this critical dimension decreases 
as more internal dimensions are considered. We learn from figure |2| that for 3 < n < 6 
the total dimension of the spacetime should be at least D = 10 to allow a continuous 
decay of an unstable black brane to the non-uniform state. This estimate is somewhat 
marginal since as we noted above in the T 1 case the actual critical dimension is 13.5 while 
the estimate is around 12.5, namely one less. This indicates that even though the estimate 
for D* is closer to 10 than in the T 1 case, the actual calculations should be performed in 
order to establish whether indeed D* = 10 or 11. 

In addition, one should exercise care using the approximate mass (|4.1| ) since it turns 
out that above certain D and n the Schwarzschild black hole "does not fit" into the torus, 
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namely its radius becomes "too large", 2po/L > 1. In such a case the equation for ^5 ( fP| ) 
is outside its domain of validity. Explicitly, one obtains 



\D-2 / 

Hence, for the D = 10 example, the Schwarzschild black hole would not fit into the torus 
for n = 5,6 being slightly larger, by 3 and 8 percents respectively, than L (see also [pS|]). 
Actually, there is a little more space in the compact dimension than appears first due to 
the "Archimedes law for caged black holes" |2^, |25j which tells us that a black hole "repels" 
an amount of space around it that is proportional to its own size. The effect vanishes in 
the large d limit, but may cause an appreciable correction for 5 < D < 10. It would be 
interesting to generalize its computation in |25| from T 1 compactification to T n . 

Another interesting estimator can be gotten by considering the Schwarzschild radius 
po for a would be black- hole with mass pgl- One finds 

D-3 _ Gd (d - 2)Q d _ 2 d -3 u on 

where ro is the Schwarzschild radius of the critical brane. For square torus Go/Gd = L n , 
and taking the large d limit where kcL = r o(2vr/L) ~ yfd one gets 

(D - 3) log (^) = j log(D) + 0(D, d), (4.4) 

from which we see that for large D (and independently of d) the black hole cannot fit in 
the compact dimension at the GL point, and therefore the black brane must decay into a 
different end-state, presumably a non-uniform black brane. This estimator can be improved 
a little by replacing the Schwarzschild radius with the radius in conformal coordinates 
Po —> Ph = po/2 2 /( d ~ 3 ^ (see [[0|), and by incorporating the "Archimedes effect", but these 
improvements would not change the large D result above. 
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A. Determination of the critical mass ijlgl{D) 

In this appendix we briefly review the first order calculations of a Lorentzian perturbation 
theory around the marginally static GL-string, B±, in D = d+ 1 dimensions. In fact, the 
first order is precisely what one needs to determine the onset of instability. An appropriate 
perturbation theory was first developed in 5D by Gubser [p]], perfected by Wiseman [11, 12] 
in 6D and generalized in ||] to an arbitrary dimension. 

The most general ansatz for static black string solutions can be written as 

ds 2 = -e 2A fdt 2 + e 2B {r l dr 2 + dz 2 ) + e 2C r 2 dn 2 



D-3' 

/ = l-(r /r) D -\ (A.l) 
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where A, B and C depend on r, z only. When these functions vanish the metric becomes 
that of a static uniform black string. We will work here in units such that the horizon 
located at ro = 1. 

Let us expand the metric functions around the marginally static uniform solution in 
powers of some small parameter 7 A such that in the limit A — ► the perturbed string joins 
the GL-point 

A = a(r)\cos(kz)+0(\ 2 ), B = b(r)X cos(kz)+0(X 2 ), C = c(r)A cos(kz)+0(X 2 ). (A.2) 

The wavenumber k is related to the asymptotic length of the compact circle by k = 2tt/L. 

Upon substituting this expansion into the Einstein equations one obtains a set of ODEs 
for a, b and c. It turns out that in our case b can obtained algebraically from 

u _ rfa + 2(D-3) fc+2rf (a' + (D - 3) d) 

b ~ 2(D-3)/ + r/' ' (A - 3) 

while for a and c there is an eigenvalue problem to solve the set of two coupled linear 
second order ODEs 

2(D-3)f + 3rf , (D — 3) f , 2 
-fa — a c+ka = 0, (A.4) 

_ „ _ 4 (D-3) f (4 -£ + (£> -3) f)-rf (4 (D - 3) f + rf) , 

r (2 (D-3) f + rf) °' 
f (2[2(4-D) + (D-3)f]+rf) 2 (D - A) f 



2(£>-3)/ + r/') r (2 (D-3) f + rf) 

c + k 2 c = 0, (A.5) 



2(D-4)f _ , u2 



r (2 (Z>-3) / + r/') 
subject to the regularity boundary conditions at the horizon, r = ro 

2 , D-3 { k 2 



a' 



k 2 a^ ( -2a + 2c 



3{D-4) 3 V D 



k 2 

c' = 2a-2c+- -c, (A.6) 

D - 4 

and at the infinity, r oo, where a and c must vanish. 

These equations do not have singular points other than horizon and the spatial infinity 



unlike the Euclidian perturbation equation (2.10). However, there is a price to pay 



instead of a single second order equation, as in a Euclidean theory, one has to deal with 
two coupled equations. 

In Q the equations (A.4- A^ ) were solved numerically for various D's. The wavenum- 



ber k obtained in this way is precisely that of the marginally tachyonic mode. (Indeed, the 
static method that is employed here converges only for this mode. For k other than kcL 
there is no perturbed static solution since above kcL the string is non static being GL- 
unstable, while below koL there is no static perturbed solution either, since perturbations 
must decay.) The calculated values of the critical wavenumber, koLi are listed in Table |TJ. 

7 A good candidate for such a parameter is the dilatonic charge or the tension along the z- direction [ju^, , 
though for our purposes A even does not have to be specified at this stage. 
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D 


5 


6 


7 


8 


9 


10 


11 


12 


kcL 


.876 


1.27 


1.58 


1.85 


2.09 


2.30 


2.50 


2.69 


D 


13 


14 


15 


16 


20 


30 


50 


100 


kcL 


2.87 


3.03 


3.19 


3.34 


3.89 


5.06 


6.72 


9.75 



Table 1: Numerically computed static mode wavenumbers kcL in units of r Q . 
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